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Large parity violating longitudinal single-spin asymmetries = — 0.86^q ^ an d -Af, = 

0.88^q are observed for inclusive high transverse momentum electrons and positrons in polar- 
ized p + p collisions at a center of mass energy of ^/s = 500 GeV with the PHENIX detector at 
RHIC. These e come mainly from the decay of W ± and Z° bosons, and their asymmetries di- 
rectly demonstrate parity violation in the couplings of the W ± to the light quarks. The observed 
electron and positron yields were used to estimate W ± boson production cross sections for the e 
channels of a{pp^W + X) x BR(W + ^>e + v £ ) = 144.1 ± 21.2(stat)^ 4 3 (syst) ± 15%{norm) pb, and 
a(pp^W-X) x BR(W~ ->e~v e ) = 31.7 ± 12.1(stat)tg° 2 1 (syst) ± l5%(norm) pb. 

PACS numbers: 14.20.Dh, 25.40.Ep, 13.85.Ni, 13.88.+e 



Determining the contributions of the partons to the 
spin of the proton is a crucial element in our under- 
standing of QCD [IHSj • Polarized inclusive deep inelastic 
scattering experiments (DIS) have measured the combi- 
nation of valence and sea quark, and gluon helicity dis- 
tributions (Ag + Aq, Ag) [31 Ej. Analysis of polarized 
semi-inclusive DIS experiments [SHI] have determined the 
individual flavor separated Aq and Ag by connecting fi- 
nal state hadrons with quark flavors using fragmentation 
functions. Collisions of longitudinally polarized protons 
at high energies allow study of Ag [U HI |9], and can 
provide complementary measurements of Au, Au, Ad, 
Ad [TD|[TT]. In particular, the production of W ± couples 
only the left-handed quarks and right-handed antiquarks 
(uLd~R — > W + and — > W~), so the asymmetry 
of the W yield from flipping the helicity of a polarized 
proton is sensitive to the flavor dependence of Ag and 
Ag. Production of the W occurs at a scale where higher- 
order QCD corrections can be evaluated reliably, and it 
is free from uncertainties in fragmentation functions by 
detecting leptons from W decay O H3] . 

The first observations of VF-boson production in po- 
larized p + p collisions, and direct demonstration of the 
parity- violating coupling of the W to the light quarks are 
reported here by PHENIX and in a companion paper by 
STAR 14 for y/s = 500 GeV at the Relativistic Heavy 
Ion Collider. The 2009 PHENIX data are from polarized 
pp — >• e + X, where the e with px > 30 GeV/c come 
mainly from W and Z decays. 

The PHENIX detector has been described in detail 
elsewhere (T3J. This analysis is based on data col- 
lected with two central arm spectrometers, each cover- 
ing \A<fi\ < ir/2 in azimuth and \rj\ < 0.35 in pseu- 
dorapidity, which surround the central axial magnetic 
field. The bend angle of charged tracks is determined 
by drift chambers outside the magnetic field starting at 
a radius of 2.02 m from the beamline. The longitudi- 
nal position z of the track is measured as it exits the 
drift chamber in a layer of pad chambers 2.46 m from 
the beamline with a spatial resolution of a z = 1.7 mm. 



The electromagnetic calorimeter, located at a radial dis- 
tance of ~5 m from the beam line, is used to measure 
the energy, position, and time of flight of electrons. In 
this analysis, the pt dependence of the reconstructed tt° 
and rj mass peaks was used to confirm the energy scale 
and linearity to within 2.5%. The px dependence of the 
peak widths was used to determine the energy resolution 
a E /E = 8.1%/VE(GeV) 5.0%. 

A trigger with a nominal 10 GeV threshold in the elec- 
tromagnetic calorimeter selected events for this analy- 
sis. This trigger was fully efficient for e with transverse 
momentum pt above 12 GeV/c. Charged tracks recon- 
structed in the drift chambers and the pad chambers 
which match the calorimeter cluster with \A<f)\ < 0.01 
were used to reconstruct the z position of the event ver- 
tex. Only events with \z\ < 30 cm that are well within 
the acceptance of the central arm spectrometers were an- 
alyzed. Loose cuts on the time of flight measured by the 
calorimeter and energy-momentum matching suppressed 
accidental matches and cosmic rays. 

The analyzed data sample corresponds to an inte- 
grated luminosity of 8.6 pb -1 , which was determined 
from beam-beam counter coincidences and corrected for a 
small (6%) effect from multiple collisions per beam cross- 
ing. The beam-beam counters are two 64 channel quartz 
Cerenkov counters ±1.44 m from the center of the de- 
tector and cover a pseudorapidity range of3.1< <3.9. 
The cross section for coincidences within \z\ < 30 cm was 
found to be 32.5 ± 3.2 mb from the van der Meer scan 
technique [15] . 

The resulting yield of positive and negative electron 
candidates is shown in Fig. [l] where pr has been deter- 
mined from the calorimeter cluster energy. The charge 
sign is determined from the bend angle, a, measured 
in the drift chamber, and the nominal transverse beam 
position. The angular resolution and stability of beam 
position were monitored by frequent runs with no mag- 
netic field. The resolution a a was typically about 1.1 
mr, to be compared to a 2.3 mr bend angle for 40 GeV/c 
tracks. The variation in the average transverse beam po- 
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FIG. 1: (color online) The spectra of positive (upper panel) 
and negative (lower panel) candidates before (solid histogram) 
and after (dashed histogram) an isolation cut. The estimated 
background bands are also shown. The computation of the 
background before the isolation cut is described in the text. 
The background band after the isolation cut is computed 
by scaling the background before the isolation cut by the 
isolation cut efficiency measured in the background region 
(12 < p T < 20 GeV/c). 



sition measured by reconstruction of the primary vertex 
in these runs was within ±300 /im, and did not affect the 
charge determination. The probability of charge misiden- 
tification at 40 GeV/c was estimated to be less than 2%. 

In addition to e ± from W and Z decay, this sam- 
ple of events contains various backgrounds. The dom- 
inant backgrounds were photon conversions before the 
drift chamber and charged hadrons. These were esti- 
mated using the raw calorimeter cluster distribution and 
the charged pion spectra predicted by perturbative QCD 
convoluted with the hadronic response of the calorime- 
ter tuned to reproduce test beam data. This calculated 
background was normalized to the measured spectrum 
in the region 12 < p T < 20 GeV/c and extrapolated to 
higher px- Electrons from heavy flavor decay were es- 
timated from an FONLL calculation [T7] which agrees 
well with the prompt electron measurement at y/s = 200 
GeV|18]. PYTHIA 19 was used to estimate the contri- 
butions of electrons with pt > 30 GeV/c from r lepton 
decays of W and Z bosons. These two components were 
found to be negligible. The background bands in Fig. [T] 
include uncertainties in the photon conversion probabil- 
ity, the background normalization, and the background 
extrapolation to px >30 GeV/c. 

The acceptance was calculated from Monte Carlo sim- 




T (GeV/cj° 



FIG. 2: (color online) Background subtracted spectra of 
positron (upper panel) and electron (lower panel) candidates 
taken from all counts compared to the spectrum of W and Z 
decays from an NLO calculation 12, 13] (see text). The gray 
bands reflect the range of background estimates. 



ulation and corrected for time dependent variations in 
tracking and trigger efficiency measured with a minimum 
bias data sample triggered by beam-beam counter coin- 
cidences. The reconstruction efficiency for tracks within 
the geometrical acceptance of the spectrometer is the 
product of the active areas of the drift chamber and the 
pad chamber (54%), a fiducial cut on the calorimeter 
to ensure the showers are not near an edge (81%), and a 
fiducial cut on the drift chamber to avoid charge misiden- 
tification (85%). The reconstruction efficiency was not 
Pt dependent for px >30 GeV/c. 

Figure [2] shows the background subtracted signal in 
the acceptance of the detector for positive and negative 
charges respectively compared to the spectrum predicted 
by NLO calculations [TU [131 HO] normalized for the in- 
tegrated luminosity, corrected for the detector efficiency 
and acceptance, and smeared by the energy resolution of 
the calorimeter. The yield measured by counting events 
in the signal (30 < pt < 50 GeV/c) region is compared 
to the predicted yield in Table [I] The measured yield in 
the signal region is consistent with the predictions of the 
NLO and NNLO [20] calculations. 

To extract the W ± production cross section, we es- 
timated the Z fraction in our sample with the NLO 
and NNLO calculations and MRST [2T] and MSTW (22] 
PDFs. The contribution from Z decays is approximately 
-7% for W+ and -30% for W~ . The fraction of the 
total cross section within \y\ < 0.35, pr > 30 GeV/c, 
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TABLE I: Comparison of measured cross sections for electrons 
and positrons with 30 < pr < 50 GeV/c from W and Z decays 
with NLO [12] and NNLO [20] calculations. The first error is 
statistical; the second error is systematic from the uncertainty 
in the background; the third error is a 15% normalization 
uncertainty due to the luminosity (10%), multiple collision 
(5%), and acceptance and efficiency uncertainties (10%). 



Lepton 


^(30<pJ<50GeV/c)| B=0 
Data NLO 


[Pb] 

NNLO 


e+ 
e~ 

e + and e~ 


50.2 ± 7.2+}>;l ± 15% 
9.7 ±3.712.5 ± 15% 
59.9 ± 8.lt!;J ± 15% 


43.2 
11.3 
54.5 


46.8 
13.5 
60.3 



and |A(/>| < ir is estimated to be ~11% of positrons 
from W + and ~7.5% of electrons from W~ from these 
calculations. The variation of the calculation is small 
compared to other sources of systematic uncertainty. 
With these corrections, a(pp —¥ W + X) x BR(W + — » 
e+v e ) = 144.1 ± 21.2(stat)l:J 4 3(syst) ± 15%(norm) pb, 
and <j(pp -> W~ X) x BR(W~ -> e~i/ e ) = 31.7 ± 
12.1(stat)tg° 2 1 (syst) ± 15% (norm) pb, where BR is the 
branching ratio. These are shown in Fig. [3] and compared 
to published Tevatron and SppS data. 

In order to determine the longitudinal spin asymmetry 
with a sample of W decays with minimal background con- 
tamination, two additional requirements were imposed on 
the candidate events. An isolation cut requiring the sum 
of cluster energies in the calorimeter and transverse mo- 
menta measured in the drift chamber be less than 2 GeV 
in a cone with a radius in 77 and <f> of 0.5 around the candi- 
date track was used to remove remaining events with jets. 
About 80% of the signal is kept, while the background 
is reduced by a factor ~ 4 as shown in Fig. [T] The sec- 
ond cut is to reject tracks with |a| < 1 mr, which reduces 
charge misidentification to negligible levels. There are 
42 candidate W + + Z a decays to positrons with a back- 
ground of 1.7 ± 1.0 and 13 candidate W~ + Z° decays 
to electrons with a background of 1.6 ± 1.0 events within 
30<pr<50 GeV/c after these two additional cuts. 

The measured asymmetry is given by 

= N+-R-N~ 
6L N++R-N- ( ' 

where N + is the number of events from a beam of positive 
helicity and N~ is the number of events from a beam of 
negative helicity, and R is the ratio of the luminosity 
for the positive and the negative helicity beams. The 
longitudinal spin asymmetry is then calculated from the 
measured asymmetry according to 

A, - ^ C) 

where P is the beam polarization and I? is a dilution 
correction to account for the remaining background in 
the signal region. 
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Theory curves: FEWZ and MSTW08 NLO PDF's^. 
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FIG. 3: (color online) Inclusive cross sections for W leptonic 
decay channel of this measurement and pp measurements [23L 
I26| . Statistical and systematic uncertainties were added here 
in quadrature. Curves represent theory calculations. 



The luminosity-weighted polarization was 0.38 ± 0.03 
and 0.40 ± 0.04 in the two beams. In RHIC, both beams 
are bunched, and the bunch helicity alternates almost 
every crossing to reduce systematic effects. The relative 
luminosities of the four helicity combinations were mea- 
sured by the beam-beam counters, and were all within 1% 
of each other. To treat the low statistics data properly, 
a likelihood function created from the four spin sorted 
yields corrected by the relative luminosity was used to 
determine the single-spin asymmetry within its physical 
range [-1,1]. 

The measured asymmetries are shown in Table [n] for 
tracks in the background (12 < pr < 20 GeV/c) and 
signal (30 < pt < 50 GeV/c) regions. For tracks in 
the background region, cl was found to be zero within 
uncertainties. Significant asymmetries were observed for 
tracks in the signal region. The dilution corrections of 
1.04 ± 0.03 and 1.14 ± 0.10 for positive and negative 
charges, respectively, were applied. 

Figure [4] compares measured longitudinal single-spin 
asymmetries to estimates based on a sample of polar- 
ized PDFs extracted from fits of DIS and semi-inclusive 
DIS data. The experimental results are consistent with 
the theoretical calculations at 6-15% confidence level for 
Al and at 20-37% for A% . The observed asymme- 
tries are sensitive to the polarized quark densities at 
x ~ Myj I ^fs ~ 0.16, and directly demonstrate the parity 
violating coupling between W bosons and light quarks. 

In summary, we presented first measurements of pro- 
duction cross section and nonzero parity violating asym- 
metry in W and Z production in polarized p+p collisions 
at y/s = 500 GeV. The results are found to be consistent 
with theoretical expectations and similar measurements 
of Af QI]. RHIC luminosity and PHENIX detector 
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TABLE II: Longitudinal single-spin asymmetries 



Sample 



A%(W + Z) 68%CL 95%CL 



Bkgrnd + -0.015 ± 0.04 
Signal + -0.31 ±0.10 



-0.86 



-1,-0.56] [-1,-0.16] 



Bkgrnd 0.025 ± 0.04 

Signal - 0.29 ± 0.20 +0.88 



[0.17, 1] [-0.60, 1] 





w* + z° 




: > 30 GeV/c 




:|yj<0.35 




* 

1 











:w"+z° 




: p* > 30 GeV/c 




:|yj<0.35 








DSSV 
DNS kkp 
DNS kre 
GRSV std 
GRSV val 





FIG. 4: (color online) Longitudinal single-spin asymmetries 
for electrons and positrons from W and Z decays. The error 
bars represent 68% CL. The theoretical curves are calculated 
using NLO with different polarized PDFs [T2"] . 



upgrades in progress will make it possible in the future 
to significantly reduce the uncertainties for and to 
extend the measurement to forward rapidity, which will 
improve our knowledge of flavor separated quark and an- 
tiquark helicity distributions. 
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